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Abstract 

A Lax system in three variables is presented, two equations of which form the 
Lax pair of the stationary Davey-Stewartson n equation. With certain nonhnear con- 
straints, the full integrability condition of this Lax system contains the hyperbolic 
Nizhnik-Novikov-Veselov equation and its standard Lax pair. The Darboux transfor- 
mation for the Davey-Stewartson n equation is used to solve the hyperbolic Nizhnik- 
Novikov-Veselov equation. Using Darboux transformation, global n-soliton solutions 
are obtained. It is proved that each n-soliton solution approaches zero uniformly and 
exponentially at spatial infinity and is asymptotic to lumps of peaks at temporal 
infinity. 

1 Introduction 

The Nizhnik-Novikov-Veselov (NNV) equation [T6l[T71[l9] is an important 2+1 dimensional 
integrable equation which is a natural generalization of the KdV equation to 2+1 dimensions. 
It is useful in both mechanics and differential geometry [HI [I2]. The NNV equation has 
been solved by various methods such as inverse scattering [2], bilinear method [18], bilinear 
Backlund transformation [8] , binary Darboux transformation [H] and so on [H El [71 El |13l |15] . 
However, one can not construct the usual Darboux transformation (without integration) 
because the principal part of the first equation of its Lax pair is two dimensional wave 
operator or Laplace operator. 

Starting from the idea of nonlinearization [3], many high dimensional integrable systems 
were reduced to lower dimensional ones so that interesting solutions like soliton solutions 
and quasi-periodic solutions can be obtained from lower dimensional systems. Especially, 
the KP equation [U [10], the DSI equation and the 2 + 1 dimensional iV-wave equation [21] 
were related to some 1+1 dimensional AKNS systems. Following this idea, in this paper, 
we present a Lax system of three variables, two equations of which form the Lax pair of the 
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stationary Davey-Stewartson n (DS n) equation. With the nonhnear constraints ([11]), the 
full integrability condition of this Lax system contains the hyperbolic NNV equation and its 
standard Lax pair. 

The DSn equation has a Darboux transformation without integration. With the relations 
given by (IHl), the Darboux transformation for DSII equation is used to solve the hyperbolic 
NNV equation. This Darboux transformation without integration is more suitable for sym- 
bolic calculation than the known binary Darboux transformation. 

It is well known that DSI equation has solutions approaching zero exponentially at spatial 
infinity, but DSn equation has not. However, we get soliton solution u of the hyperbolic NNV 
equation from that of the stationary DSII equation so that u approaches zero exponentially 
at spatial infinity. This is possible because the solution u of the hyperbolic NNV equation 
is given by i{g — g) as in (fT4|) . not /, the solution of the stationary DSII equation. These 
soliton solutions are different from the known one derived by binary Darboux transformation 
or bilinear method etc. and the behavior of the solutions is more complicated. 

In Section [21 after reviewing the hyperbolic NNV equation and the stationary DSII equa- 
tion together with their standard Lax pairs, a new Lax system f[TO|) is presented in which an 
extra equation is added to the standard Lax pair of the stationary DSII equation. With the 
nonlinear constraints ([HI), the integrability condition of this Lax system includes both the 
hyperbolic NNV equation and its standard Lax pair. The Darboux transformation for the 
new Lax system is given in Section [3] and the general expression of multi-soliton solutions is 
presented in Section [H In Section [5l the explicit expressions and behavior of single-soliton 
solutions are discussed. In Section [6l it is proved that each n-soliton solution approaches 
zero uniformly and exponentially at spatial infinity. In Section [TJ it is proved that each 
n-soliton solution is asymptotic to 'n? lumps of peaks at temporal infinity. Finally, some 
linear algebraic lemmas are presented in the Appendix. 

2 Hyperbolic Nizhnik-Novikov-Veselov equation and 
Davey-Stewartson IE equation 

The hyperbolic NNV equation is 

Ut = U^^^ + Ur^r,r, + 'i{uv)^ + d>{uw)^ri, 



which has a Lax pair 



ft = /ccc + fmv + + ^^fv 



By taking the new coordinates x = ^ — rj,y = ^ + rj, the hyperbolic NNV equation ([T|) 
becomes 

Ut = 2Uyyy + QU^^y + ?,{u{V + W)) y + ?,{u{V - W))^, 

[dy - d^)v = {dy + d^)u, {dy + d^)w = {dy - d^)u, 
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and the Lax pair ([2]) becomes 



fyy fxx 0, 

ft = 2fyyy + GU^^y + 3(f + W)fy + 3 {V " W) f ^ 



(4) 



On the other hand, the DSn equation is 



which has a Lax pair 



-i/r = fxx - fyy " ^9 " 9)f, 

(dy-id,)g = id,-idym'), 

^y = iJ-^^ + P-^, 



(5) 



(6) 



where 



-1 / V -/ / V - i/j/ g 

If (/,(?) is independent of r, ([5]) becomes the stationary DSn equation 

fxx - fyy - Kg - g)f = 0, 

{dy-id.)g = {d.-idy){\f\'). 

Taking ?/,r) = $(2;, ?/)e^'^^'^ in Qj, we get the Lax pair for ([8]) as 

% = U^^ + P^, 

2iA^$ = 2iJ$^.^. + 2P$^. + Q$ 



(9) 



The first equation of (jl]) and the first equation of ([8]) are similar, and the second equation 
of (jl]) is of order 3. Hence we introduce an extra equation to the Lax pair iQ so that the 
whole system becomes 

^y = M(a)$ = iJ$^ + P$, 

2iA2$ = L{d)^ = 2iM^^ + 2P^^ + Q$, (10) 
= N{d)^ = 1617$,,, + 16P$,, + + 

where J, P, Q are given by ([7]), 
i? = 4 



5 = 2 



3^? + i|/r 4/,-2i/, 
-4/,-2i/, 3^-i|/|2 

3g., + 2i/7. + //, - //, Qfxx - 2if,y - i{g - g)f + 2|/|V 
-6/... - 2i^ + i{g - g)f-2\f\'f 3g, - 2ifU + //, - ffy 



and L{d), M{d) and N{d) refer to differential operators with respect to x whose coefficients 
are 2x2 matrices, d = d^- 

The integrabihty conditions of ( fTOl) include the following equations: 



fyy fxx + uf 0, 

ft = '2fyyy + Q f ^^y + + W) f y + ?,{V " W) f ^ , 



(12) 



|2n 



{dy-id^)g={d,-ldy){\fY 

'2^9^ '^9xxx ~l~ "^f fxxy ~l~ '^f fxxy ~l~ '^^f fxxx ~l~ "^^f fxxx (13) 

+2(/. - i/,)/., + 2(/, - i/,)/,, + 2(i/, + 2fy)f,, + 2(i/, + 2/,)/,, 
+ (2|/p - i(^7 - ^))(|/P), + (6i|/p + (^7 - ^))(|/|\. - 2\f\% + 6igg^, 



where 

w = i((7-^), v = 2\f\' + ig + -g), w = 2\f\' - (g + g). (14) 

Note that f|T2|) is exactly the same as the original Lax pair (jl]) of the hyperbolic NNV 
equation. By direct calculation, we know that {u, v, w) satisfies the hyperbolic NNV equation 
([3]) provided that / and g satisfy ( fT2l) -( fT4l) . Therefore, explicit solutions of the hyperbolic 
NNV equation can be obtained from those of (|T2l)-f|T^. 

Clearly, the solutions of (fT2l)-( |T^ are only part of those of the hyperbolic NNV equation. 
However, they include some interesting ones which will be shown in the rest of this paper. 

3 Darboux transformation 

The binary Darboux transformation for the hyperbolic NNV equation is well-known [Tl]. 
Integrations are needed in constructing explicit solutions. However, for DSn equation, usual 
Darboux transformation without integration is known. This Darboux transformation is 
simpler than the binary Darboux transformation for the hyperbolic NNV equation, and can 
be easily used to the stationary DSn equation so that explicit solutions of the hyperbolic 
NNV equation can be constructed. 

Note that the coefficients of L{d), M{d), N{d) satisfy 



where 



U,P,Q,R,Se^ (15) 



S = {A is a 2 X 2 matrix | KAK~^ = A} = { \ _ \ a,b E C} , (16) 




K= (^^ ^ j. That is, L((9), M{d), N{d) satisfy 

KL{d)K~^ = L{d), KM{d)K'^ = M{d), KN{d)K-^ = N{d). (17) 
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Hence, if $ = f ^ j is a solution of ffTOj) with A = Aq, then = i j is a solution of 



([TO]) with A = ±iAo. 

The Darboux transformation of arbitrary order is constructed as follows [SJ [20] . Suppose 

G{d) = d" + G^{x, y, t)d^-' + ■■■ + y, t) (18) 

is a Darboux operator for ([TOl) . i.e., there exist L'{d), M'{d), N'{d) which have the same 
form as L{d), M{d), N{d) with / and g replaced by certain /' and g', such that $' = G{d)^ 
satisfies 

W = L'{d)^', ^'y = M'{d)<^', <l>; = iV'(9)$'. (19) 

If so, G{d) satisfies 

L'{d)G{d) = G{d)L{d), 

M'{d)G{d) = G{d)M{d) + Gy{d), (20) 
N'{d)G{d) = G{d)N{d) + Gtid). 

Since L{d), M{d), N{d) satisfy the relations ([IT]), and L'{d), M'{d), N'{d) satisfy the similar 
relations 

KL'{d)K'^ = L\d), KM'{d)K~^ = M\d), KN'{d)K-^ = N'{d), (21) 
we want that G{d) satisfies KG{d)K'^ = G{d). Write 

Denote L'{d) = 2iJd'^ + 2P'd + Q', then the first equation of ([20]) leads to 



{2iJd^ + 2P'd + g')(^" + GiS"-! + ■ ■ ■ + G„) 
= (9" + + ■ ■ ■ + Gn)i2Ud^ + 2Pd + g), 

in which the coefficients of d'"'~^^ and 9" give 

P' = P-i[J,Gi], 

Q' = Q- 2i[J, G2] - 2[P, Gi] + 2i[J, Gi]Gi + 2nP, - 4iJGi,. 
Hence, after the action of Darboux transformation. 



(23) 



(24) 



/' 


= f- 2i6i, 








9' 


= g - 4iai_^ - 


2(/6i - 


fbi)- 


-4i|6i|2, 


u' 


= M + 8|6ip - 


4i(/6i - 


-fhi) 


+ 4(ai + ai)a;. 


v' 


= i; + 8|6i|2 - 


4i(/6i - 


-fhi) 


- 4i(ai - ai)^. 


w' 


= w + 8|6ip - 


-4i(/6i 


-fbi 


) +4i(ai - ai)^ 



(25) 
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Now take n distinct complex numbers Ai, ■ ■ ■ , A„ with \j = fij + iuj {fi/s and uj^s are real). 
Let = 



be a column solution of ( fTOD with A = \j, then ^n+i = K^i — 



■ n+J 



IS 



Vj J _ V 

a solution of (fTOj) with A = ±iAj {j = 1 ■ ■ ■ ,n). The Darboux transformation is determined 
by the system of linear algebraic equations 



G{d)^,=0 (j = l,---,2n) 



if it has a unique solution 
Denote 













( ^ 












[ Vn ) 


, a = 


\ an ) 


, b = 


[ hn ) 



then (l26l) becomes 



where 




A = ( 5"-^^ ■ ■ ■ ^ ) , 5 = ( d'^-^ri 
has a unique solution if and only if det T 7^ 0. 



(26) 



(27) 

(28) 

(29) 
(30) 



4 Expression of soliton solutions 

For zero seed solution u = v = w = f = g = 0, (ITU]) becomes 

= AV$, = iJ$,, $i = 1617$. 
with $ = (^, 77)-^. Hence take 



(1) , ■ (1) (1) • (1) tO\ (2) , . (2) (2) . (2) 

(J-J/ P' +10-. I —p. — icr \ „ (^)/ p +10-. I —p. — i(T ' 

J — Kj 3 + e ^3 3 ), rjj = Kj [e'^i ^ j + e '^j ^ ^ 



where 



p'P = Re(AjX + iXjH 

Pf = Re (iA 

- ^ + iXjy + 16iA|t) + ajo^ = 



-r i/vjy + 16iA|t) + = yUj-a; - Uj-y + 16(z/| - 3fi'jUj)t + pfj, 
yjX + XjU - 16A^t) + Pjo^ = -VjX + /ij?/ - - 3yUjZ/|)t + p^o\ 



a)^' = Im(A 
.(2) 



af ^ = Im(iAjX + X^y - IQX^t) + 



-UjX + fijy — 16(pi- — 3fiji^j)t 
= UjX + jj.jy + 16(/i^ - 3jj,ji'j)t + (T^^^ 



fijX + Vjy + lQ{v^ 



(31) 



(32) 



(33) 
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nf^ are non-zero constants, p]•o^ p}o^ ^}o^ ^jo'' real constants. By solving a/s and 
6j's from fl28l) . the Darboux transformation (1251) gives the n-soliton solution 



Jl) J2) ^(1) ^(2) 



/ = -2i6i, 5( = -4iai,^ - 4i|6ip, 



M = 8|6ip + 4(ai + ai; 



- 4i(ai - ai)^, w = 8|6ip + 4i(ai - ai)^. 



(34) 



Hereafter we omit the primes on /, u, v, w for those obtained by the action of Darboux 
transformation. 



Let Kr, 



. Denote 



-7] 



for j > k, then 
Let 

n = 



^ — ( Rn-l -O KnRn-1- ) • 



d^-\ d^-\ Rn^3-0 KnRn-1-0 



(35) 
(36) 

(37) 



Theorem 1 When det T 7^ 0, the multi-soliton solution u of the hyperbolic NNV equation 
given by (34^ can be written as 

^ Re det H 
u = —i 



Proof. Solved from (|28|) by Cramer rule, 



ai = -(detT)-i Rn-2-0 K^R 



n-l---0 



6l = -(detT)-^ i?„_1...0 d^C KnRn^2-, 



ai + ai 

= -(detT)-i( 
= -(detT)-i( 



d'^C Rn-2-0 KnRn-l-.Q 
d"'C Rn-2-0 KnRn-l-.Q 



<9"C Rn-2-0 KnRn-l-0 
Rn-1-0 Knd^C, KnRn-2-0 



-(detT)-i(detT), = -tr(T-^Tj, 

+ di,, = - tr(r-ir,.,.) + tr((r-ir,)2). 



(38) 

(39) 
(40) 

(41) 
(42) 
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Denote K 



[ai,---,a„)'^, h 



Let h 



Ikxk 
0(n— fc)xA; J 

= (61, ■ ■ ■ then 




be the solution of Th 



d"'^^( where a 



o-i^x + 0'i,x — ~ tr 



+ tr 



—a 




6 


-b 


-b 


—a 


—a 


In-1 b 





-b 


-a 


In-1 



-a 
2 





In-2 



(43) 



+ + ai + ai — 02 — 02 — 2|6i| 



According to ([3] 
Let c?= -(detT)2(a2 + 



02 J 



f^"C Rn-2-0 KnRn-1- 



u = 8Re(a^ + Si — 02). 
then by Cramer rule, 
2 







+ detT 



71-3--0 



+ detT 

KnRn-l--0 



n-2---0 



n-fLn-l-.-O 



(44) 



(45) 



Using Laplace expansion of det LI, d 



det n. Hence 

Red „ Re det U 



u = —i 



{det ry 



(det ry 



(46) 



The theorem is proved. 



Remark 1 According to LemmalE o/Hl detT > holds everywhere. However, detT > 
may not hold everywhere when the parameters pf^ and crjg^ take some special values, as will 
shown in the next section for single soliton solution. On the other hand, det T > holds 
everywhere in generic case, which will he shown here. 

The Darboux operator G{d) of order n can he constructed hy composing n Darhoux op- 
erators of order one as follows. For given Ai, ■ ■ ■ , A„ and as above, let Hj = 
($j, $„,+-,) j{= l,---,n). If det Hi ^ 0, then Ai((9) = d — Hi^Hi^ is a Darhoux oper- 
ator of order one. It transforms {u,v,w,f,g) to (u^^\v^^\w^^\ f^^\ g^^^) and transforms 
Hj to H^ = Ai{d)Hj = Hj,x - Hi,xHi^Hj {j = 2,3,---,n). Agam, if det H^^^ ^ 0, 
then A2id) = d- H^^liH^y^ is a Darhoux operator of order one for the Lax pair with 
(u«,t;«,ti;«, /«,(?«) . It transforms (u«, /«, ^7^) to (m^^), ^(2)^ ^(2)^(2)^ ^{2)) 
and transforms Hf^ to H^ = A2{d)Hf^ = HfJ, - H^^iiH?)'^ Hf^ (j = 3, 4, ■ ■ ■ , n). Con- 



tinuing this process, we get H^ 
According to f2^. 



(fc) 



{k 



j = k + 1, 



n) and Aj {d) {j 



G(a)=A„(a)A„_i(9)---Ai(a), 
detT = det{Hi^-^^) det{Ht~i^^) 



■det{Hi^^)det{Hi] 



(47) 
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Hence det T ^ if all det h!-^ ^ V 0. 




Suppose H^/-''> = \^ I , then det iff = 1^]'"'^^ + = z/ and 

only if ^j''^^^ = and T]f ^^^ = /ioW simultaneously. For fixed j , this gives a system of four 
real equations 

Re^j^-i) = 0, Imej'"'^ = 0, Reri^^-'^ = 0, Imr/j^"'^ = (48) 

for three real variables x,y,t. It has no solution unless the parameters p^'^J and ajo'' (j = 
1, ■ ■ ■ , n; k = 1,2) take special values. This shows that det T > holds everywhere for generic 
pf^ and <y^jQ . Therefore, the multi-soliton solution u is global for generic p^g^ and cy^jQ . 



5 Single soliton solution 

By taking n = 1, the single soliton can be obtained as 



where 



(50) 



B = {H^l^ f cosh(2pf )) + (/s:f V cosh(2pf )) + (/cS'V cos(2ai'^) + (/«f ^ cos(2frf )) 
A = -{pi - z/2)(4'^)4cosh(2p^')) cos(2ai'^) - 2piz/i(4'))4 sinh(2pS')) sm{2a[^'>) 
-{pi + z/2)(/cS^')2(42))2 sinh(2pf ^) sin(2(Tf ^) 

+ {pl - z/2)(4'Vcosh(2pf^) cos(2ai'^) + 2piz/i(/€p))^sinh(2pf^) sin(2af)) 
+ {pl + z/2)(/tS^))2(42))2gij^i^(2pf)) sin(2ai'^) 

+{i^i-ui){{n?yr-{n^:yn 

The solution is singular if B = 0, i.e. + = 0. This is equivalent to p^^^ = p^^^ = 0, 
2a\^^ = jn + 7r/2, 2a^^ = kn + it /2 for certain integers j and fc. In contrast, the solution is 
global if and only if + |?7ip 7^ everywhere, i.e. the parameters satisfy 

/^i(pSo^ - +kn + n/2) + u,{p'S + - jn - n/2) ^ 0. (51) 

We always suppose (15T|) is satisfied, which is equivalent to detT 7^ 0. 

When pI uf, the solution u approaches zero exponentially at spatial infinity, and the 
peaks appear when neither p^i^ nor p^i^ is large. Hence the center of the lump of peaks 
locates near pi^^ = and p^^ = 0, i.e., 

(1) I (2) (1) I (2) 

X = QApiUit 2 2 ' y = 16(Pi + Vl)t 2 2 ■ (^2) 

Pi — Ui Pi — Ui 
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Figure 1: Single soliton solution u: Ai = 2 + 0.5i, Hi'^ = 1, n^i^ = 1.2, p^Q = 0, p^iQ 
4'^ = cr^^ = 0, t = 1. 



The solutions are shown in Figure [T] and Figure [2] for different parameters. The figure of 
the solution contains a lump of peaks rather than a single peak, and the shape depends on 

the angle arctan — between the straight lines p^^"* = and pf^ = 0. Nevertheless, we 

2/iiz/i 

still call it single soliton solution because it is generated from the zero solution by Darboux 

transformation, and the peaks in the solution never separate. 

Note that although u is localized, v and w are not. 

1 / \ 1 1 / rmT\ 

If 1^1 = Pi 0, the solution is invariant when ix,y] is changed to x H ,y H 

^ 2/ii 2pi^ 

for any integer m. Hence the solution is periodic. Moreover, p^i^ + pf^ = p^iQ + p^g^. The 
peaks appear when neither pi^^ nor p^i^ is large. Hence the peaks lie near the straight line 

(1) _ (2) 

X — y ~ 32p?t H = 0. The solution is shown in Figure El 

2pi 

Similarly, the solution is also periodic if ui = — pi ^ 0. 

6 Localization of the solutions 

In this section, we will prove that the multi-soliton solutions approach zero uniformly and 
exponentially at spatial infinity. In order to get global solutions, we always suppose det T 7^ 
everywhere, which is true for generic parameters p^g'* and ajo^ (j = 1, • • • , n; = 1, 2). 

Note that the solution of fl28l) is invariant if both and r]j (for fixed j) are multiplied 
by a common function. Let 

[ vj if m < m 

T = diag(a;i, ■ ■ ■ , w^, cui, ■ ■ ■ , u;„). (54) 



10 



Figure 2: Single soliton solution u: Ai = 1.1 + 0.9i, Ki'^ — 1, k^i^ — 1.2, p^^ — 0, p^Q 
= = 0, i = 1. 



Figure 3: Periodic solution u: Ai = 1 + i, ac^^'* = 1, k^'^ — 1.2, p^^ — 0, p^Q — 1, ctiq = 
ai? = 0, i = 1. 
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Let T = T ^T, then the norm of each entry of T cannot exceed 1. Although T is not 
continuous, |detT| is continuous. 

Let X = rcosO, y = rsin6'. Since pf^^s and pf^^s depend on x and y hnearly, we can 
write, for /c = 1, 2, 

pf (r cos e, r sin 6, t) = ef {e)af {e)r + (3^^ , (55) 

where e^J'\6) = ±1 (j = 1,- ■■,n) so that a^''\6) > 0. Here the variable t is omitted in 
af and ef \e). 



Clearly a] s are continuous functions. Note also that Ej (6) is not well-defined when 



af\e) = 0. 

Theorem 2 Suppose Ai, ■ ■ • , An are distinct non-zero complex numbers such that Xj ^ ±iAz 
for all j,l = 1, ■ ■ ■ ,n. u is the n-soliton solution given by ^EE)- Then for fixed t, there are 
positive constants tq, x one? C such that 

|u(r COS0, r sin 0, t)| < Ce"^^ (56) 

for r > ro and all e'^ G S^. Hence u{r cos 6, r sin6,t) —>■ uniformly and exponentially as 
r +00. 

Proof. The proof is divided into four steps. 

Step 1: Obtain the asymptotic behavior of ^/s and rj/s. 

Let Xj = pj + Wj where /i/s and zz/s are real, then pj ^ ±.Uj for all j = 1, ■ ■ ■ , n. 

Let Z(^) = {e''' G \ tan^ = e} fore = ±1, Z = Z^^^^UZ^-^\ If e'^ G Z^'\ then by dSS]), 

af{e) = af\e) = -^\pj - euj\ > 0, ef\9) = eef{d) and ef{d){pj - euj) cos^ > for all 

j = 1, ■ ■ ■ ,n. If e'^ G 5'^\Z, then af\6) = \pj cos 6 — uj sin 6^1, af\0) = \ —J^j cos 6 + pj sm6\ 

with af\e) 7^ af\9). 
For S G (0,7r/4), define 



there exists e'^" G Z(^) such that 161 - 9o\ < S] (e = ±1), 

l^-^ol > 5/2 for all e'^" G Z}. ^ ^ 



af\9) > u, af\e) > uj, ef\e) = eef\9), ef\9){pj - ev^) cos^ > if e'^ G , 

(0) 



Then U fij U df^ = S^, and there exists 6 G (0, 7r/4) and cu > such that 

«j'^(0) >cu, «f (0) >cu, ef\e) 

\a'j'\e) - af\e)\ > 00 ife'^en 

(58) 

For e'^« G Z^^) and e'^ G l^i^^ with 16* - 6*01 < 5, ef V) is a constant, 
«f^(0)-af (6^) 

= ef'^ (6) {pj cos - Uj sin 0) - ef {e){-iyj cos ^ + pj sin ^) (59) 
= eef\6){pj + ez/j) cos 6^(6 — tan6'). 
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Hence, if af\e) > af\e) for e'^ G fif \^^'^ with < - < 5, then af\e) < 
e'^ G Qg^\z^'^^ with — 5 < 6^ — < 0, and vise versa. 
Recall that 



a 



(2) 



Vj 



' (e)r+e(.i) (e)^^ (0,r) 



J2), 



2q('^ {e)r-2ef^ {9)13^^ -2;sf {e)af^ (9,r) 
^(2)g"f' {e)r+ef {e)pf + ief {e)af (e,r) ^ ^-2af {e)r-2sf (e)/3f > -218^ {9)af (e,r)^ 



^(2) (^E^'(e)XjX+E^'mjy-^^'iS)>^t ^ ^-ief\e)X,x-ef\e)X,y+l&ef\e)x:p^ 



Let 



When r — > +cxd, the following limits hold uniformly. 
For e'^ G with aj^^) > V), 



-(1) 



0. 



For e'^ G with < ^^'^(e), 



(1), 



,,(2), 



For e'^ G with a'^-^O) > af>{9) 



„(2), 



„(2), 



For e'^ G with aj-^^) < (^Y'i^) 



Step 2: There exists ro > and cq > such that det T > cq when r > tq. 
When e'^ G Z(^) (e = ±1), a^p{e) = af\9) for j = 1, ■ ■ ■ , n. dSHD implies 



\mr)\ 



as r ^ +00. By ((621), 

det f(^,r) 



n 

|7,{e)|>i 



A{d,r) B{9,r) 
,r) A{e,r) 



where y4(6',r) and B{6,r) are n x n matrices, whose entries are 



+ 0(1), 
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and 0(1) refers to the terms which tend to zero as r — > +00. Let 



r= diag(ri(^,r),---,r„(e,r)). 



(69) 



By Lemma [3l of [Xf there exist ri > and Ci > such that detT > ci for e'^ G Z^^"^ and 
r > ri. 

When e'^ G = ±1) with < ^ - ^0 < ^ (^0 e ef\e) = eef{e). Suppose 



af{e) > a)"' (6) for j 



(2), 



m and q;j-^''(^^) < <yY'{6) for j = m + 1, ■ ■ ■ , n. By ( l62l) and 



„(2) 



det T(^,r) 



C{e,r) D{e) 
-B{e,r) A{e) 

-D{e) c{e,r) 



+0(1) 



A{e) B{e,r) 

D{0) -C{e,r) 

-B{e,r) A{e) 

c{e,r) D{e) 



+ 0(1) 



(70) 



where A{6) and B{6,r) are m x n matrices, C{6,r) and D{6) are {n — m) x n matrices, 
whose entries are given by 



m 



A,,ie) = isf\e)x,)\ B,,ie,r) = iie4\e)x,)%ie,r) (j = i, 

Dj,{e) = {-ief\9)\,)\ -C,,{9,r) = {eef\9)Xj)\-Y,{9,rr') (j = m + 1, 



(71) 



(72) 



and 0(1) refers to the terms which tend to zero uniformly as r — > +00. 
Let 

A = (4'^(e)Ai, ■ ■ ■ ,£«(e)A„, -ieSi(^^) Wi, ■ ■ • , -i4'H^)An), 
r = diag(ri(^^, r), ■ ■ ■ , r), -F™+i(0, r)-\ ■ ■ • , -F„(^, r)-i). 

By LemmaHofEl there exist > and C2 > such that det f > C2 for e'^ G 1]^+^^ ur]^"^^\Z 
with < — 6*0 < 5 and r > r2. 

Similarly, when e'^ G 0^+^^ U with -5 < ^ - ^0 < (^0 G ^), there exist rg > 

and C3 > such that det T > C3 for r > rs. 



When e'^ G suppose aj.'^(0) > aY^{9) for j = l,---,m and aj'^(^) < a)"' (9) for 



„(2) 



„(2)/ 



j = m + 1, ■ ■ ■ , n, then, by (IMl) and (!65|) . 



lim det T'(6', r) 

r— >+oo 



A(^) 







A{9) 








D(0) 




m 








A(^) 







m 


-D{9) 










m 



A{9) 

m 



(73) 



holds uniformly, where A{9) is an m x n matrix, D{9) is an {n — m)x n matrix, whose entries 
are given by 



A,,{9) = {ef {9)\^)\ (j = l,---,m), 
Djk{9) = [-\ef{9)X,f (j = m + 1, ■ ■ ■ , n). 



(74) 
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Using the condition Xj ^ ±iA; and the property of Vandermonde determinant, we know that 
there exist > and C4 > such that detT > C4 for all e'^ G and r > r^. 

Let ro = max(ri, r2, r^, r4), cq = min(ci, C2, C3, C4), then for any e'^ G 5*^, det T > cq when 
r > tq. 

f f 

Step 3: Denote H = o n, then lim Re det 11 = for any fixed e'^ G S^. 

When e'^ G Z^^\ considering ([62]), ([631) and £f^(0) = eeJ-^V), let A = diag(Ai, A2) with 

Ai= diag(eS'H^)Ai,---,e«(^)A„), 
A2 = diag(-i££S'^(e)Ai, ■ ■ ■ , -i££«(^)A„) 

C = (1, ■ ■ ■ , 1, -er'^i, ■ ■ ■ , -C'^n)"^, (76) 



then A2 = ieAi, and H— ]^ l7j(^)l) ^ as r ^ +00 where 7j(^)'s are defined by 

|7,W|>1 

( 1^ and is defined by (I122p . According to Lemma HI of Rl Re det 11^ = 0, which leads 
to lim Re det fl = 0. 

r— >+oo 

When e'^ G S^\Z, suppose af\9) > af {6) for j = 1, ■ ■ ■ ,m and af\e) < af\9) for 
j =m + l,---,n. By (leij) and (l65|) . 

^0, e7i(9'=r/,-(i£«(^^)A,)S) -0 (j = l,---,m), 

^-^ld%-{ef\9)\,)%)^0, r^7i(5% -(ief (^)A,)S) -0 (j = m + 1, ■ ■ ■ , n) 

(77) 

as r — > +00 since ^r]j — > for j = 1, ■ ■ ■ , m and r]j — for j = m + 1, ■ ■ ■ , n. 
Let A = diag(Ai,A2) with 

Ai = diag(4'^(^)Ai, ■ ■ ■,Eg\e)Xm.eS:l^mm+i, ■ ■ ■ ,4'n^)A„), 
A2 = diag(-i£S'H^)Ai, ■ ■ ■ , -ie(^\e)Xm, -i£2+i(^) Wi, • • ■ , -i4'H^)An), 

C = (l,---,l,0,---,0,0,---,0,-l,---,-l)^ (79) 



(7J 



m n—m m n—m 

then A2 = iAi, and 11 11^ as r — > +00 where ll'^ is defined by (11221) . According to 
Lemma H] in |Al we have lim Re det 11 = Re det ll'^ = 0. 

r— >+oo 



Till now, we have proved that detT has a uniform positive lower bound for all 9, and 

lim — — — — = for any fixed 6. 
^-+00 (detT)2 ^ 

Step 4: lim -— = uniformly for all e'^ G S*^ as r — +00. 

^ r^+cx) (detT) 2 ^ 

Note that — — ——r is of form '^f ' \ where 
(detT)^ g{0,r) 

f{e,r) = ^e"^(^)'-+^^(^), g{e,r) = Y^ef'^'^'>+'^^'^ (80) 
j=i j=i 
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are real- valued functions of {0,r), C(j{6), (3j{6), 7j(^), Sj{6) are (complex valued) continuous 

~ ~ ~ /(^,^) 

functions of 9. Let a(9) = max Re a, (6'), b(9) = max Redj(9). Since lim , = 

i<i<mi i<i<m2 ' ■'^ ^ r^+oo g(9,r) 

for any fixed ^, the real continuous function a{9) — b{9) < achieves its maximum — x < 
on the compact set S^. We have known that detT has a uniform positive lower bound as 
r > ro (ro is independent of 9), so has 



"12 



{/3,(e)-6(e))r+5,(e) 



(81) 



Hence 



mo 



g{9,r) 



< e 



mi 



{a,(e)-a(e))r+7,{e) 



< Ce 



-Xr 



E< 



,(/3,(e)-6{0))r+5,(e) 



(82) 



as r > ro where C is a constant independent of 9. The theorem is proved. 



7 Asymptotic behavior of the solutions as t ^ oo 

In this section, the asymptotic behavior of the n-soliton solutions as t ^ oo will be discussed. 
In order to do so, we consider the problem in a moving frame. Let x = xo + 9it, y = yo + 
where {9i, 92) is the velocity of the moving frame, and (xq, |/o) is the coordinate in the moving 
frame with this velocity. Then 

pf{x^ + ^it, 1/0 + e2t, t) = efaft + (5f\ 
pf\xo + 9,t, yo + 92t, t) = efaft + (3f 

where ef'^ = ±1 (j = 1, ■ ■ ■ , n; k = 1,2) so that af^ > 0. Write Aj = fij + ivj {j = 1, ■ ■ ■ ,n) 



(84) 



where /x/s and Uj's are real, then according to fl33|) . 

= jj^Oi - uj92 + 16(z/3 - 3/i|z/,), 
ef'af^ = -i^j9i + fij92 - 16(;U^ - 3/ijZ/|), 

and 

(1)^,(1) _ .(2)^,(2) _ , ^1R/„2 „2 



(/ij + z/j)(^i -92 + 16(/ij - + z//)), 

e^^a^^ + ej'^ttf = (/x, - i/,)(^i + ^2 - 16(/i| + AfijU, + u])). 



jl) Jl) , J2) (2) 



(85) 



Theorem 3 Suppose Xj = fij + iz/j 7^ (j = 1, ■ ■ ■ ,n) are distinct complex numbers where 
fij 's and Uj 's are real numbers, such that 

jij 7^ ±z/j for all j, 

p] + 4:pjVj + z/2 ^ p1 + 4fiiui + vf for all j ^l, 
p] - AfijUj + ^ pf - 4:HiUi + vf for all j ^1. 
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u is the n-soliton solution given by (E 
6^2^, t) = except when 



Then for bounded {xo,yo), lim u{xo + 9it,yo + 



t—*CX) 



02 = Sifif + Afiiui + uf + fi'j - 4:iJ,jiyj + z/2 
(j,Z = l,2,---,n). 



01 = 8(/if + Af^iui + uf - ii] + AiijVj - z/|), 



Therefore, as t oo, u has at most n x n lumps of peaks which move in the above velocities 
{61,62) respectively. 



Proof. 



We will always suppose that (01, ^2) does not satisfy ([87]). Then, by (ESD, ^ 
— af\ Moreover, we only consider the limit t +00. The conclusion is the 



whenever 



same for t —00. 

The proof is divided into three steps. 

Step 1: Obtain the asymptotic behavior of ^j's and ?7j's. 



Suppose af^ > af^ for j = 1, ■ ■ ■ , m; af^ < aY' for j = m + 1, 



j = P + 1, 



> of 
. n. Then 



(1) 



af ^ for 



0, i-\d\ 
0, ^f(d\-{\ef\,)\)- 



-^0 (j = l,---,m), 
(j = m + 1, ■ ■ ■ 
(j = p+ 1, . . . ,n) 



(88) 



as r ^ +00 where Si, ■ ■ 
TjJ^^j ^ for j = m + 1, 



Sp are any constants, since ^- ^rjj — > for j 



,P- 



Now we prove that p can only take n — 1 or n. If p < n 



2, then ef^ef = 

(1)^(2) 



must hold for any j 7^ / with p + 1 < j,l < n since both sides equal ±1. If e)j Sj 



3 



-(1)^,(1) ^(2)^,(2) 



a 



then e) 'a 

contradicts condition ( l86l) . If ^f^^f^ 



(1)^,(1) 



hold simultaneously where e 



-e 



(1)^(2) 



then e 



(2)^(2) 



(1)^,(1) 



eel " 



hold simultaneously where e = ef^ef\ This contradicts the assumption that 



, m and 



efef. This 
= 0, e?a? + 



'3 



?i, 62) does not satisfy fl87|) . Hence only p = norp = n — lis possible. 
Step 2: There exists to > and c > such that detT > c for t > to- 
When p = n — 1, 



\in{6.r)\ 



max(|en(0,r)|,|r/„(0,r)|) 



7o 



I (1)1 e'^^fl^'^ 



/, (1), .(l)fl(l) 

max(|Kn |e^" . 



(2)Uel^)/3(^) 



(89) 



Denote 



n-1 



1 \ 
1 / 



(i-j+l)xn 



(90) 
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for j < I, then, for p = n — 1, (188|) leads to 



det T = 7^ 







C'r/„r(i42)A„) 










T// • (2) \ • (2) \ 

V(4'^Ai,---,e«A„) 
V{ie[^^i\m+i-, ■ ■ ■ 








ie„_lA„_iJ 



as t — s> +00. Let 

A = diag(£:i"'"'*Ai, ■ ■ ■ , e^-'A™, — iel^^^Am+i, ■ ■ ■ , 

r= diag(o,---,o,cV), ^ = 4M'\ 



— le 



X P-(i)x ^ 



(91) 
(92) 



then we get hm inf det T > by Lemma 151 of [XI 

t— >+oo 



Similarly, when p = n, 



lim detT 

t— >+oo 



r(eWAi,---,e«A„) 













\/(4'^Ai,---,£«A^) 



i^m+lAm+l; ■ ■ ■ ; i^i, ^A^ 



(93) 



det \/(£S'^Ai, A 



Using the condition A,- ^ ±iA;, we get lim inf detT > since the Vandermonde determinant 



IS non-zero. 



Step 3: Denote LI 



T 



T 



n. If the velocity (6*1, 6*2) does not satisfy (187|) . then 



lim Re det LI = 0. 

From dHH]), let A = diag(Ai,A2) with 



A^= diag(4^)Ai, A 

A2 = diag(-i£:f^SiAi, ■ ■ ■ , -ie^^^s^A^, -ie^+iAm+i, ■ ■ ■ , -i4^^Ap, 
— iffp^iAp+i, ■ ■ ■ , — ie^^^An), 



(94) 



Figure 4: 3 x 3 soliton solution u: Ai = 2 + 0.5 i, A2 = 2.5 + 0.4i, A3 



,(2) 



1.2, p« = 0, p);' =0,a 



(2) 



.(1) 



0, 



(j = 1,2,3), t = l. 



3 + 0.3i, nf = 1, 



C = (1, ■ ■ ■ , 1, 0, ■ ■ ■ , 0, 1, ■ ■ ■ , 1, 0, ■ ■ ■ , 0, -1, ■ ■ ■ , -1, -^p+iVp+i, ^VnV, 



(95) 



p—m 



n—p 



p—m 



then n — for p 



n 



1 and n — 11^ for p = n as t 



by 

for j 



-00 where 70 is defined 

(1)^(2) 



and n is defined by (11221) . A2 = i^Ai holds for e = ±1 if and only if sj = ee, e 



p, and ef^ = eey for j = p+1, ■ ■ ■ ,n. Since sj {j = 1, ■ ■ ■ ,p) can be arbitrary, 

r{l)p(2) 



-(1) 



e can be taken as ±1 arbitrarily, and p = n or n — 1, we have A2 = ieAi by taking e = e^n'^n 
for p = n — 1, and either £ = lor£: = — 1 for p = n. 

According to Lemma HI of Rl Re det 11"^ = 0. Hence lim Redetll = if (01,62) does 

r— ►+00 

not satisfy (I57|) . The theorem is proved. 



Remark 2 1/(61,62) satisfies ( [^^ , then a^''' = af \ a[^^ = a] ' , ' = ey , ef — — 

j 7^ / with p + 1 < j,l < n. Hence there is no common e = ±1 such that ef^ = ee\^'^ holds 
for alii = p + 1, ■ ■ ■ ,n, which contradicts the condition A2 = ieAi in Lemma^ oflAl In fact, 
the solution does not tend to zero in this case, which can be seen in the following example. 

As an example, a 3 x 3 soliton is shown in Figure HI in which there are 9 lumps of peaks. 
The local behavior of each lump of peaks is still complicated and one of which is shown in 
Figure [51 



(2) (2) 



-(1) ^(2) 



(1) / 

—el for 



A Some linear algebraic lemmas 

Lemma 1 Suppose X and Y are 2n x r and 2n x (2n — r) matrices respectively, then 

X K„Y = Y K„X 



(96) 
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Figure 5: Local behavior of one lump of peaks in the 3x3 soliton solution u. 



where = ( ° 



Proof. 

X KnY )k-' = {KnX -Y)(^^^ o" ) = ( ^ ) ■ ^^^^ 

The lemma is obtained by taking the determinants on both sides. 

Lemma 2 Suppose T = \ _ _ | where A and B are n x n matrices, then det T > 0. 

\ -B A J 

Proof. By Lemma [T], det T is real. First suppose both A and B are invertible, then 

detT = det A det{A + BA-^B) = \ det A\'^ det{I + A^A-^B). (98) 

Let = A~^B. Suppose A is an eigenvalue of NN, v G R^" is a vector in the corresponding 
root space TZx, i.e. {NN — XI)"^v = for certain positive integer m. Then 



{NN - XI)"'{Nv) = N{NN - A/)™w = 0. (99) 

Hence Nv G TZx- If A is a non-real eigenvalue of NN of multiplicity k, the multiplicity of A 
is also k. 

Now suppose A < is an eigenvalue of NN, Tlx = Vi © ■ ■ ■ ©Kn where Vi, - ■ ■ ,Vm are irre- 
ducible invariant subspaces. Suppose Vi = span{C, {NN — XI)C, ■ ■ ■ , {NN — A/)™~^C} with 
{NN-XI)'^C = 0. The n C ^ lmage{NN - XI), and {NN-XI)'^NC = N{NN - A/)™C = 0, 
{NN - A/)™"iiVC = N{NN - A/)'^-iC since det ^ 0. We will prove that 

C, {NN - A/)C, ■ ■ ■ , {NN - XI)"'-\, iVC, {NN - XI)NC, {NN - XI)"''^NC (100) 



20 



are linearly independent. Suppose 



j=l j=l 



(101) 



where ai, ■ ■ ■ , a^, Pi, - ■ ■ , Pm are complex numbers. Acting (A^A^ — A/)"* ^ on both sides of 
dm]), we get 

(102) 



Then 



{NN - A/)'"-^(aiC + PiNQ = 0. 



= -aipi{NN - A/)™-^iVC - X\PiWNN - XI)'^-\ 
= -aipiN{NN - A/)"^-iC - X\piWNN - XI)""^C 
= piNpiiNN - XI)'^-WC - X\Pi\'^{NN - A/)™-iC 



(103) 



= |/?i|2A^A^(ArA^ - XI)"'-^C - X\pi\^{NN - XI)'^'^C 
= \pi\\NN - XI)"'C = 0. 

Since A < and {NN — XI)"^^^( ^ 0, we have ai = Pi = 0. Continuing this process 
by acting (A^A^ - A/)""-^ ■ ■ (A^A^ - A/)° on both sides of ffTOTD respectively, we get 
ai = ■ ■ ■ = an = Pi = ■ ■ ■ = Pn = 0- This proves the linear independence of the vectors in 
dm . Let Vi = span{A^C, (NN - XI)NC ■ ■ ■ ^{NN - XI)"'-^NC}- If NC = (NN - XI)C' G 
Image(A^A^ - XI), then C = {NN - XI)N~'^(' G Image(A^A^ - XI), which contradicts the 
choice of (. Hence N( ^ Image(A^A^ — XI). Moreover, Vi is invariant and irreducible under 
the action of A^A^ — XI. Hence it must be one of Vj with 2 < j < m, which means that m 
is even (m = 2k) and Tlx = {Wi ® Wi) ® ■ ■ ■ ® {Wk © Wk) where {Wi, ■ ■ ■ ,Wk,Wi, ■ ■ ■ , Wk) 
is a permutation of (Vi, ■ ■ ■ , V2k)- Therefore, the multiplicity of each negative eigenvalue of 
A"A^ must be even. 

Thus, if the eigenvalues of A'A^ are Ai, ■ ■ ■ , X2n (multiple eigenvalues are listed repeatedly), 
then 

det T = I det det(/ + A^A^) 



detA|2( n |l + A,f n + 



1/2 



n (1 + A,) > 0. 

A,>0 



(104) 



If y4 or i? is not invertible, it is a limit of invertible matrices, and the conclusion is also true. 
The lemma is proved. 



Lemma 3 Suppose A = diag(Ai, 
such that A J 7^ iiA^ for all j,l = 1, 



V 



( XT' Ar' 

\n— 1 \n—2 



■ , Xn) where Ai 
■,n, 

1\ 



Xn are distinct complex numbers 



n,T= diag(7i,72, ■ ■ ■ ,7n). Denote 



, i?= diag((i£)"-\(i£; 



n-2 



(105) 



1/ 
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where e — ±1. Let 

( V WE 

T=\ _ 

\^ -VVE V 

then there is a positive number C depending on A only, such that det T > C . 
Proof. Denote 



r-(m) _ \ \ 

then 



Ji<---<ifc 



n n 
k=l '=1 

Since the entries of V are Vjk = A^~^, we have 

s=l 

Hence 

det T = dGtiV) deiiy + VVEV-^VVE) 
= |det(l^)|2det(/ + fl^^l^-iri^^l^-i). 

For j = 1, ■ ■ ■ , n, let Aj = e"^/'^5 where 5 = 1 if £ = 1 and 5 = i if £ = — ] 
distinct and fj ± r; 7^ for all j, / = 1, • • • , n. 



3=1 



" Ti + r. 



where 



7=1 - Ts 



A = {ajSjk) , B = {bjSjk) , M = {Mjk) , 

n n 

S=l a=l 

s^k 



Since VEV-^VEV-^ = I, we have AMB-^ = {AMB-^)-\ Hence 

det T = I det(y)|^ det(7 + TAMB-^TBM-'^A-^) 
= I det(y)p(detM)-i det(M + fMr) 
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since F, A, B are diagonal matrices. 

Suppose Re(ri), • • • , YLe{rm) < 0, Re(r^+i), • • • , Re(r„) > since fj + rj for all j. 
Write 

1 \ / Mil Mi2 



M + TMr 



(115) 

where Mn and A^n, M12 and A^i2, M22 and A^22 are m x m, m x (n — m) , (n — m) x {n — m) 
matrices respectively. Then Mn and Nu are negative definite Hermitian matrices, so are 
Mil ^11 I -M22 and iV22 arc positive definite Hermitian matrices. 
Let 

Ti = diag(7i, • • • , 7m), ^2^ diag(7m+i, • • • , 7n), (116) 



then 



Hence 



det(-iVii) = det(-Mii + rt(-Mii)ri) > det(-Mn), 
det(Ar22) = det(M22 + r^MaaLa) > det(M22). 

(-1)™ det(M + f Mr) = det(-iVii) det(Ar22 + Nl,{-Nii)-^Ni2) 
> IdetMiil • |detM22|. 



(117) 



(118) 



On the other hand, 



(-l)™det M = det(-Mii) det(M22 + M*2(-Mii)-iMi2) 
> IdetMiil • |detM22| > 0, 



(119) 



hence 

, _ , , ,„JdetMii| ■ |detM22| 

det T > det V ^^- , ' ' , (120 

' |detM| ^ ^ 

The lemma is proved. 

Lemma 4 Let A = diag(Ai, A2) where Ai and A2 are n x n diagonal matrices satisfying 
A2 = isAi where e — ±1. Let ( be a2n dimensional column vector. For j > k, denote 

Rf,.., = (A^C A^-\,...,A'C). (121) 

Let 

nA^i' ^"C A"-^C A'^-'C Rts-o K^Rti-o 

where = ( ^ " ) , then det H^ is purely imaginary. 
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Proof. A2 = ieAi is equivalent to AKn = —ieKnA. Denote = det A^. Multiplying 
row j of by —\j and adding it to row 2n + j {j = 1, ■ ■ ■ , 2n), we get 



d' 



AX A"-^C A"-'C Rts-o K^Rn-i. 















R'n-2-1 ^^^nRn-.-O R'n-2-0 ^nRn-1-0 



by using AKn = —ieKnA. Adding columns 2n + 2, - ■ ■ ,3n — 1 to columns 4, 
multiplying columns 3n+ 1, ■ ■ ■ , 4n — 1 by —ie and adding them to columns n + 3, 
we get 

C..0 KnA^-'C KnRt2-o 
ieir„A"C 

By moving the columns, 



d' 



Rn-2-Q ^nRn~l-0 



d' 



pA pA 
J^n-Q -"■n-'l'n-2-0 












pA jy- pA 

^n-n'n-2-0 




Rt 



KnA^~\ 



K 



(123) 

,n + 1, 
2n+l, 

(124) 



(125) 



which is purely imaginary according to Lemma [H The lemma is proved. 
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